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ABSTRACT 

The effects of heat transfer on the kinetic energy tendencies of a 
perturbation superimposed on a baSic zonal current are investigated. The 
rate at which kinetic energy is gained in both the adiabatic and diabatic 
models increases with time for medium wave lengths, but the rate of in- 
crease is less for the diabatic case. At the longer wave lengths the 
energy change for the diabatic model shows a decrease with time becoming 
negative toward the end of the time interval investigated. The adiabatic 
model, in contrast, has zero energy tendency at long wave lengths. The 
rate of kinetic energy change increases as larger values of thermal wind 
are selected, and for the adiabatic model the wave length range for which 
energy changes occur also expands with increasing thermal wind. The 
energy gain initially is a maximum when the thermal wave lags the pressure 
wave by a phase angle of 90 degrees and the maximum tendency to lose 


energy occurs when the lag is 270 degrees. 
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the relative vorticity 
the wind velocity vector 


the average value of the Coriolis parameter and equal to the 
the value at latitude 45N 


the gradient operator in pressure coordinates 
the acceleration of gravity 

the height of the pressure surface 

the temperature 


the total derivative of pressure with respect to time (vertical 


velocity) 
a WE es 
the static stability parameter (55 Cp 


the Rossby parameter 
the gas constant for dry air 
the Jacobian operator 

mt 2 
a measure of static stability ( psR ) 
the specific heat for constant pressure 


the stream function (24 (2+22) ) 


the thickness function ( 4 (2,-2,)) until specified 
as 7 = 3.14159 be 


the heat gain or loss 

the atmospheric pressure 
the wave number ( a 

the wave length 

the phase speed of the waves 
time 


the basic zonal wind 


= the basic thermal wind 


TABLE OF SYMBOLS 
(Continued) 


the relative phase angle between temperature and pressure wave 
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I. Introduction 

To date simple baroclinic models have been unable to improve signif- 
icantly on the results obtained by one parameter models in the prognosis 
of mid-tropospheric circulation. Recent investigations have shown that 
the incorporation of previously neglected terms, which were felt to be 
negligible for short range forecasting, baci ee pailt deren: Berea Marae the 
characteristics of synoptic waves produced by baroclinic forecasting 
models, thus bringing these waves into closer proximity with snyoptic 
waves observed in nature. 

‘In one of these investigations, Haltiner [2 ] considered the influ- 
ence of sensible heat exchange on the properties of harmonic perturba- 
tions superimposed on a zonal current. The study revealed that this 
type of heating reduces the instability of short and medium waves, but 
that the instability of long waves is increased. 

The purpose of this work is to compare the kinetic energy changes 
averaged over the interval of one wave length in a baroclinic model with 
heat transfer present to one when heating is absent. In particular, the 
manner in which heating alters the instability of long waves will be 
examined. A simple two-level quasi-geostrophic model is constructed « 
and an equation expressing the time rate of change of kinetic energy 
in a harmonic perturbation superimposed on a zonal current is derived. 
The kinetic energy iS computed numerically for the range of wave lengths 
corresponding to synoptically observed waves. At each wave length com- 
putations are made for varying values of thermal wind and phase lag be- 
tween the pressure wave and the temperature wave. The change of kinetic 
energy with time is evaluated in time steps of three hours over a thirty 


hour period and the value of the tendency for the heating and the non- 


heating model are compared. 

In addition, computations are made of the values of C;/L, which is 
proportional to the growth rate of the wave, and the values of phase 
speed to verify that the results obtained from the simple model used are 
in qualitative agreement with those obtained by Haltiner. 

Interest throughout the study is confined to unstable waves only 


and solutions which lead to neutral waves are eliminated. 


Il. Mathematical Development 

The atmosphere is represented by a four-layer baroclinic model 
using the 750 and 250 millibar levels as data levels and centering atten- 
tion at 500 millibars. Since the effects of latitude variation are well 
known, latitude is held fixed at 45 N. The vertical velocity is assumed 
to be zero at both the top and the bottom of the model atmosphere, the 
effects of friction at the earth's surface are neglected, and the static 
stability of the atmosphere is taken to be, constant. 

The vorticity equation is applied at the levels corresponding to 
250 and 750 millibars in the following form: 


a) v-V(I+7)= PV. (1) 


The continuity equation is used to Sct for the divergence in the 
above equation and then the derivative of the vertical velocity with 
respect to pressure is approximated by a finite difference ratio. Geo- 
strophic approximations are substituted for the vorticity and wind. The 


sum and difference of this equation at the data levels are computed with 


the following results: 


raha ( 2). T(z (Ue) me 4, Ve. “)/- E 


(2) 
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At the middle level (500 mb) the First Law of Thermodynamics is 


applied in the following form: 


arly V-VT +35, 5 > : +: (4) 


SBE a 


Jt 


where Dar _ 
3) a ( op Ce Zz o 


The temperature at this level is approximated using the equation of 


State and the hydrostatic equation to give 
R (2,- 2.) 


This approximation is used in (4) and used to replace ® in (3). 


With the following substitutions: 
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2+ # (2 3) TT = 7, (2-8) Mo 4 fe 


fo RS, 

fic, Weaken Sgenieaecoele 
(2) +T(pv*y) + 3 (1, V*r) +p 3? | ra 
v7 (3) A(R): 3 (4, Vr) + 317, vp) 
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(5) 


Let the basic flow and the surface be considered to be in thermal 
equilibrium. Similar to Dods [1] , the exchange of sensible heat is 
then assumed to be proportional to the difference between the surface 
perturbation air temperature (TD and the temperature of the ground (T_) 


and a decreasing power of pressure in the following manner: 


te . (7, EY (7) 
where Gis a “hie constant. 
Since ridges in the thermal wave will correspond to heat sinks for 


the air, and troughs will be heat sources, let the transfer of heat be 


proportional and opposite in sign to the perturbation thickness. There- 
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fore 


ik 
Km’ =~). at (8) 


where K iS a proportionality constant with units sec7l, 
To linearize the system (5) and (6), assume that the solutions are 


of the form 


tk (eet iR(y-e 
y- bs ¥-ct ) qué, i (y + | 


be PG) HOS) we Ty) + w (Ke) 


Dropping terms involving products of perturbation quantities leads 


to a homogeneous linear system MN=0, where M is the matrix 


“k(e-U)+ LRA ~ik? 4 


Lk Ula A?) ik (et och ak 


(9) 


and N is a column matrix of db and Wf. In order to have non-trivial 
solutions to this system of equations the determinant of the coefficient 
matrix M must vanish. Setting this determinant equal to zero yields the 
frequency equation 
5 2 
)+(; 2 ~ 2i6kK-b°k*)) 
rek™ a) 4b kus Rk?) J (10) 
2 2\ \V/ 
_ wk )} . 
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— . aif t= _atr 
d4 a4 p= 9 ay . 
It can be seen from this equation that the presence of heating in- 


sures that the phase velocity will be a complex number regardless of 


dl 


the wave length. If the heating is taken to be zero, it can be shown 
that the phase velocity will be complex only for a certain range of wave 
lengths. The width of this range will depend on the Strength of the 
thermal wind and the value of the static stability. 

Using the primary roots of the complex number under the radical sign, 
the phase velocity has two solutions: C, when using the positive sign in 
front of the radical, and C_ when using the negative sign. The complete 


solutions for ip and 77 are therefore given by 


p y* ik (x-¢,t) cklx-et) 
=e 2 


me po 2 (11) 
+ ckle-e,t4) ck (v-ce) 
Y= 1% & + TT, & (12) 


and 


where the primes are omitted and henceforth yp and 7f will be perturba- 
tion quantities unless otherwise specified. 


Next take as initial conditions at time t=0 


p = ey T+ Be 


where A and B are real numbers and ™ is a relative phase difference 


‘(hire ) 


which is positive when the thermal field lags the stream field. From 
(11) and (12) it follows immediately that at t=0 
t = = CO ( is) 
Y), a y) = 4A 7, t TT, = be 
Along with (13), use of the linearized form of (5), once for C, and 
; : R a = We = : 
once for C_, yields four equations in “ Y * Tt, » and W, - Letting 


Cp equal the Rossby wave speed, the following solutions result: 


M74 
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The above solutions for y," > Y~ » a and T, are complex func- 





tions. Therefore, they can be separated into their real and imaginary 


parts 


i = Ine re bt 


~ r = 
and similarly for " , W772... ana *. 
o Oo r) 
As mentioned previously, the phase velocity can be a complex num- 
ber depending on whether heating is present and the values of the thermal 


wind and static stability. The phase velocity can then be represented by 


= aecbelaeed) DY a es 
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The kinetic energy of the perturbation may be expressed in terms of 
the square of the velocity of the perturbation wind. Since yp is inde- 
pendent of y and considering a unit mass, the kinetic energy per unit 


2 o) 
KE b= 2% 


A solution for V is obtained by substituting into (11) the solutions 
for C, and C_ as given in (15) and then differentiating (11) with respect 
to X. The real part of V is squared, multiplied by one half (1/2) and 
integrated over an interval of one wave length to give the average 


kinetic energy of the perturbation over one wave length as follows: 


: I" Vide = BL BOR (gate yas) «LO (dee) 


hk | CUE Ue oar ¥ Yor cima) | 


: (17) 
wk [MCAS Wie cinats UE tor oaar') | 


where ¥= kat and 7— in this and subsequent equations is Wf = 3.14159. 
The time rate of change of this kinetic energy is then obtained by 


differentiating (17) with respect to time with the following result :: 


2(¥ dy = aia b (e+ Yer) +24 (Hoe + or ) | 
k atbOR a (doe + z)- ~Qd (Yor + Vor )] 
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This energy change equation when evaluated at some time t=t, repre- 
sents the tendency of the perturbation to gain or lose kinetic energy. 
This gain or loss then shows whether the wave is amplifying or dampening. 
Since only unstable waves are studied, a solution leading to zero corres- 
ponds to a neutral wave. 

Haltiner [3] has shown that the mean energy changes over one wave- 
length for adiabatic conditions can be expressed as 


a (Ke) ro - { o.(y-) de 


Ate 
and (19) 


J ( Pe) + uy) # d 
- = fod w) oY 


where KE and PE are the mean kinetic and potential energies. The quantity 





(¢ - ¢, ) is equal to 7{ ; so after substituting the values for CG, and eC. 
into (12) and solving for / from (4), the energy transformation func- 
tions on the right hand side (19) may be determined. The kinetic energy 
changes that have been examined in this paper are thus a consequence of 


this energy transformation from potential to kinetic energy. 


15 


III. Results 

It is obvious from equation (18) that a determination of precisely 
how heating affects the stability characteristics of harmonic waves would 
be difficult. Examination of the frequency equation (10) does show that 
with heating present an instability term (term b) is always present, and 
also that heating leads to an instability term evolving from the radical 
part of the equation (term d). Term b always leads to a dampening effect 
while term d may have either dampening or amplifying characteristics. 
Therefore, the combined effects of these two terms will determine the 
energy change characteristics of the heating model. In the adiabatic 
model the energy change is dependent on only term d, which takes on a 
different value in general when heating is present. 

To arrive at quantitative results for energy changes in the two sit- 
uations, representative values of the different parameters involved are 
selected and equation (18) is solved numerically. The latitude is fixed 
at 45 degrees north latitude, the basic flow zonal wind remains constant 
at 20 m/sec, and the static stability also remains constant with "be = 
-6K/100mb and so M = 2.4 x 107!© om-2, The initial amplitude of the 
stream function and the thickness function are taken to be equal with a 


9 10 em*/sec which corresponds to a z-amplitude of 20km. 


value of 20 x l 
Since both the thermal wind and phase lag are of significant importance 
in the energy change, four values of thermal wind (5, 10, 15, and 20 
m/sec.) and 12 phase lags, 30 degrees apart, are used. 


It remains only to select a representative value for the heating 


proportionality constant. The value: for G in (7) is chosen to be 


G= lOO ERGS =" OK! een 
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and the value for r is r=2. This is similar to values used by Doos L | 
and corresponds to a total heat transfer of 63 cal per square centimeter 
per day to the entire column of air for every one degree temperature dif- 
ference between the ground and the air at the surface. To determine a 
value for K, the initial perturbation thickness is translated into a 
temperature perturbation by use of the hypsometric equation. Since the 
surface air and the ground are assumed to be initially in thermal equil- 
librium, this temperature perturbation also equals the maximum difference 
between the air and the ground. A perturbation thickness of 20 x 1019 om? / 
sec which is the selected initial value translates into a temperature of 
6.54K and corresponds to a heat transfer of 654 ergs g7l sec, A value 


for K can now be calculated from (8). The value for K at 500 millibars 


is then arrived at by using r=2 and turns out to be 


K=2269 x10” see” 

Using the above values, the first items determined were the stabil- 
ity characteristics of the waves. As mentioned in the introduction this 
is the value of the complex part of the phase velocity multiplied by k, 
where a The results are summarized in figures 1 through 4. Figure | 
1 is for no heating and a thermal wind of 10m/sec. Figure 2 is again for 
no heating and now the thermal wind is 20m/sec. Figures 3 and 4 are the 
corresponding situations for the heating model. 

To gain sufficient resolution between the values for each model at 
interested wave lengths the graphs were scaled such that the values at 
short wave lengths often exceeded the bounds of the graph. Where the 
bounds have been exceeded such that the graph plot is a straight line 
parallel to the x-axis, the actual values at these wave lengths are 


greater than those depicted on the graph. In cases where the values 


7 


are significant to the problem, special mention will be made of the values. 
It should also be mentioned that the data points are obtained at intervals 
of 1000 kilometers and between data points the computer simply makes 
straight line extrapolations. 

Examination of figures 1 through 4 reveals the following items. Heat- 
ing leads to a smaller amplification factor than no heating and a larger 
dampening factor. The values at 5000km, which are off the graph for a 
thermal wind of 20m/sec, also show that the no heating leads to a larger 
amplification and less dampening. The effect of thermal wind is readily 
seen to increase the instability characteristics even though the value of 
the b term in the heating model is independent of thermal wind. Also it 
is important to note that in the adiabatic case the instability is con- 
fined to a finite wave length interval which depends on the thermal wind 
and static stability. In contrast the presence of heating leads to an 
instability factor at all wave lengths, however the amplification ap- 
proaches zero with increasing wave length. These results are in agree- 
ment with similar computations made by Haltiner. 

Figure 5 depicts the relationship between the speed of wave propaga- 
tion in the two models. Ignoring the value at wave length zero for the 
heating case, heating leads to a decrease of phase velocity of the wave 
at all wave lengths. This difference is, however, quite small at 4-6000 
km. This also is in general agreement with results obtained by Haltiner. 

Comput ations of the time rate of change of kinetic energy show some 
interesting features. First of all, if the energy change tendency at 
time t=0 is evaluated for several different values of phase lag, the max- 
imum positive effect occurs at a phase lag of 90 degrees and the maximum 


negative change is at a phase lag of 270 degrees. The pattern is sym- 
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metrical about these two points for both models with the exception that 
for the adiabatic case the values are identically zero at phase lags of 
O and 180 degrees while in the case of heating they are small but not 
identically zero. A plot of the kinetic energy tendency at time t=0 for 
several different phase angles is shown in figures 6 and 7. As the time 
was increased from zero the phase lag of maximum positive and maximum 
negative energy change tended to increase, so that by 12 hours several 
situations (wave length and thermal wind values) showed a maximum posi-~ 
tive value at 120 degrees of lag and maximum negative at 300 degrees. 
Also, in all situations the magnitude of the positive rate of change was 
greater than that of the negative change. Again the graphs indicate that 
the presence of heating leads to kinetic energy changes at all wave 
lengths while its absence restricts the energy change to a finite inter- 
val. This phenomena remains true throughout all situations investigated. 

Another point that should be brought out at this time is that in 
the derivation of equation (18), both wave components of (11) and (12) 
were used. The energy change equation, therefore, represents the com- 
bination of the damped and amplifying components. Hence, the values of 
wave length, thermal wind, and phase lag which result in positive energy 
change represent situations when amplification will take place and vice- 
versa. 

When computations for the kinetic energy tendency for the heating 
and non-heating case were compared at t=0, the values were identical 
for corresponding situations (with the exception of phase lag O and 180 
degrees which was previously mentioned). This seemed quite surprising 
at first since equation (17) for the diabatic and adiabatic cases is 


quite different. Upon further reflection it is apparent that since 


Pe 


heating has had no time to exert an influence, the values would be the 
same. Figures 8 and 9 show the tendency at t=0 for several thermal winds. 
On the curves for a thermal wind of 20m/sec, the values at wave lengths 
5000, 6000 and 7000 kilometers can be seen to be equal in both cases. 

Next the tendency is evaluated at t=6 hours. [In both cases the 
values increased in magnitude but now an examination of figures 10 and 
11 show that the values for the heating case are smaller than those for 
no heating. When time is allowed to increase further, the difference in 
tendency values for the two cases becomes larger. Figures 12 through 15 
are for t=12 and t=24 hours and show how the difference increases. By 
24 hours the values at 5000 kilometers are off the graphs but the raw 
data still gives a higher value for the adiabatic situation. 

Values of energy change were actually tabulated out to 33 hours at 
a phase lag of 90 degrees. In the adiabatic model the tendency value 
increased out to 33 hours in all cases. When heating was present the 
value also increased with time for medium wave lengths with the exception 
that after 30 hours for a wave length of 5000 kilometers and a thermal 
wind of 5 and 10m/sec, the value started to decrease. This suggests that 
not only does the tendency increase at a slower rate when heating is in- 
troduced, but also that it reaches its maximum value Boone: At the 
longer wave lengths the value of the tendency decreases with time such 
that by 24 hours kinetic energy is actually being lost by the waves with 
wave lengths larger than 10,000 to 11,000 kilometers. It was mentioned 
earlier that the phase lag of maximum energy change increases with time 
also. This change of phase lag does not, however, alter the results that 
heating leads to smaller values of tendency as the time increases. 


So far we have examined primarily the spectrum of medium length waves 
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to determine the effects of heat transfer on the energy changes. At the 
long wave lengths the addition of heat transfer to the model leads to 
only a small rate of energy change. This, however, is significant since 
in the adiabatic case the energy change is identically zero at long wave 
lengths. That heating insured the presence of an instability term at 
all wave lengths was observed in equation (10) and again in the results 
of figures 3 and 4. These did not, however, indicate whether the effect 
was to amplify or dampen the wave. Examination of the kinetic energy 
tendency due to heating clearly shows that the effect depends on the time 
and the phase lag. For example, the series of graphs on energy change 
show that out to at least 24 hours, if the phase lag is 90 degrees, the 
heating leads to an amplifying wave at long wave lengths. Beyond 30 
hours the heating leads to a dampening effect. 

The energy tendency at the very short wave lengths was not examined 
comprehensively and no explanation is offered for the very large values 
of energy change observed. It can be seen from the graphs, however, that 
besides displaying very large values, the kinetic energy tendency may 
change abruptly from a negative to a poSitive value or vice versa as the 


time or wave length varies. 
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